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An efficient quantum algorithm for the one-dimensional Burgers equation 
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We analyze one-dimensional classical and quantum microscopic lattice-gas models governed by a 
lattice Boltzmann equation at the mesoscopic scale, achieved by ensemble averaging over microscopic 
realizations. The models are governed by the Burgers equation at the macroscopic scale, achieved 
by taking the limit where the grid size and time step both approach zero and by performing a 
perturbative Chapman-Enskog expansion. The quantum algorithm exploiting superposition and 
entanglement is more efficient than the classical one because the quantum algorithm requires less 
memory. Furthermore, its viscosity can be made arbitrarily small. 

PACS numbers: 03.67.Lx 



Here we present the simplest example where a quantum 
computer is demonstrably more efficient at numerically 
predicting the time-dependent solutions of an important 
nonlinear one-dimensional partial differential equation, 
the classical Burgers equation 



dtu + ud x u = vdlu, 



(1) 



a simplified model of turbulence and shock formation 
with flow field u(x, t), sound speed c s , and kinematic vis- 
cosity v. We shall derive ([!]) as the general effective field 
theory describing the large-scale behavior of microscopic 
one-dimensional lattice gas models with two particles per 
site. The occupation probabilites for the two particles at 
position x and at time t are denoted by p+(x,t) and 
P—(x, t), respectively. The mesoscopic kinetic transport 
dynamics is governed by the lattice Boltzmann equation 



p±(x±5x,t + 5t) =p±(x ) t)±n(p+,p-), (2) 

where denotes the collision term, a nonlinear function 
of p + and p_. The particular functional form of de- 
pends on the model type, either classical or quantum 
mechanical. For any ft, the model (J2J) always conserves 
particle number density defined as p = p+ + p-. The 
propagation speed of particles is the ratio of the lattice 
cell size to the time step interval, c = 4|. In these mod- 
els, the macroscopic flow field is 



u(x, t) = c (p(x, t) — 1) 



(3) 



The total classical computational complexity is the 
product of the lattice size N x , the ensemble of size N, 
and local resources of size g 



C A = NN x g. 



(4) 



The quantity g is the amount of resources needed to en- 
code ft for the local microscopic collisions. Since we are 
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modeling the time evolution of a classical system, the 
value of u(x, t) must be known everywhere at every time 
step whence the need for continual measurement in any 
quantum model. Therefore, the memory load factor for 
the lattice is still N x and not log 2 N x in the quantum 
case [Q. Furthermore, because of wave function collapse 
by Von Neuman projective measurement, ensemble aver- 
aging is required just as in the classical case, so a lower 
bound for the quantum computational complexity C qu is 



C q „ > NN X log 2 g. 



(5) 



For some quantum computing technologies, such as spa- 
tial nuclear magnetic resonance quantum computing 
|g, [3j, N may be counted as order unity overhead. Relat- 
ing g to the minimum number of bits per node needed to 
encode @), as a parsimonious demonstration of ([|) and 
(||) , it is known (|l|) can be modeled classically with 3 bits 
per node whereas we prove quantum mechanically only 
2 qubits per node is required. 

Multiplying ([!]) by u and integrating over all space, 
with periodic boundaries, gives a relation for energy con- 
servation where the time-rate of change of the turbulent 
kinetic energy density 9t(^) is balanced by the viscous 

dissipation e = v (ff ) 2 ~ where £ is the character- 
istic scale of the largest feature in the flow field and uc 
is the standard deviation of the turbulent kinetic energy 
or eddy velocity at that scale. 

The flow velocity, the kinematic viscosity, and the vis- 
cous dissipation quantities have the dimensions: [u] = 4, 

i 

[v] = ip-, and [e] = 5%. The dissipation scale A = (jrj 
is the smallest spatial scale where macroscopic effec- 
tive field theory ([!]) is physically applicable and the 

smallest physical velocity at the dissipation scale is the 

1 

dissipation-scale velocity u\ = (ve) 4 = j. 

The Reynolds number characterizing the fluid's non- 



to 



linearity is Re = -x«7 = 1^1 anc ^ us i n g e = ~t 

4 / 2 \ 2 4 
eliminate £, we also have Re = ^ = ( ^ J = (j) 3 - 

The computational complexity of numerically modeling 



2 



fluidic behavior at the macroscopic scale can be expressed 
as a function of Re. First, the number of grid points N x 
to sufficiently resolve the flow field down to the dissi- 

r 1 

pation scale is N x = y = Re 4 . Second, we get N by 
physically limiting |5j?±|. Since, the occupation proba- 
bilities are measured by ensemble averaging over N in- 
dependent microscopic realizations, \Sp±\ ~ ^7=, ^ ue ^° 
either classically stochastic shot noise or quantum me- 
chanically stochastic projective measurement. Following 
Orszag and Yakhot Q, the value of the statistical fluc- 
tuation Su{x 1 1) of the numerical flow field must be much 
much less than U\(x, t) of the macroscopic effective field 
theory (Q) . Using Q , we have u\ c(Sp + + <5p_ 

Hence, the ensemble size N ^> 4r ~ ^jj. 

Mach number is M = Uc/c. 

It is convenient to treat the occupation probabilities 
as a two-component field 



where a z 



1 

-] 



^ J . Then using ( |l2| ) we have 



\6 P ) 



J+ - J. 



-a z Sxd x \d) +0(e 2 



(14) 



Taking the difference of the respective components gives 



Sp + — 6p- 



J + - J. 



-Sxd x p + 0{e 2 ). 



(15) 



Similarly from (^), we obtain the second order equation: 



where the 



Std t \d) + a z 5xd x (\d) + \6p)) 
+ ^ 2 |d) + 0( £ 3 ) = 



n 
-n 



(16) 



We now take the sum of the respective components: 
Std t p + 5xd x (d + — <i_ + Sp + — 5p_) 



\p) 



(6) 



+ ^d 2 xP + O(e 3 ) = 0. 



(17) 



We expand \p) about equilibrium its value denoted \d) 
so that \p) = \d) + \5p) + 0(e 2 ), where e ~ is the 
Knudsen number. The equilibrium condition = 

leads to a tractable polynomial equation for d±, whence 
the linearized finite-difference equation 



Inserting ( |15|) into the above equation gives the gen- 
eral effective field theory for any one-dimensional two- 
particle-per-site lattice gas conserving particle number 



\p(x ±5x,t + St)) - \p{x, t)) = J\5p(x, t)) , 
where the Jacobian of the collision term is 



J : 




p—d 



j+ J- 

-j+ -J- 



(7) 



(8) 



d t p + cd x (d H 

Sx^ 
+ 2St 



-d_) 
2 



6x 2 d x (J+-J-) 



d x p 



J+ - J_ 



St ( J+ - J- 
+ l)dlp + O(e 3 ) = 0. (18) 



The left and right eigenvectors of J are 



(6 



J+ - J- 



(£l| = (l 1) 
(J + J_) 



16) 
16) 



J- - J+ 



J- 
-Ji 



1 

-1 



(9) 
(10) 



The Briegcr-Bonomi model of the Burgers equation has 
one bit per site with a 3-bit stencil where the center bit 
is updated by a stochastic Masters equation ||. This 
nonlocal classical model is equivalent to a 3-bit per site 
local lattice-gas model. We briefly discuss the classical 
Boghosian-Levermore lattice-gas model || || before pro- 
ceeding to our quantum model. The classical collision 
term is 



^ = \{P- -P-i 



2 (P++P- - 2 P+P-) 



(19) 



with associated eigenvalues Ai = and A2 = J+ — J- 



J may be rewritten as 



j = A 2 ie 2 )(6i = 



•A 



-J- 



(11) 



where a is occupation probability of an additional ran- 
dom bit used to bias the collision ||] . Taking the equilib- 
rium to be of the form d± = ^ ± A, then Q\ p=d = leads 
to a quadratic equation in A. We take the negative root 
A = a| (l — and hence 



J is singular. Nevertheless, its generalized inverse is 



M -J+ — J- 



J+ 
-J4 



-J_ 



(12) 



Now we invoke the continuum limit where Sx — > 
and St — > so |p) is a continuous and diffcrcntiable two- 
component field. We obtain a first order equation by 
Taylor expanding (Q) in x and t and keeping terms only 
terms first order in e: 



With J± = 

J+ - J- 



dp± 



d+-d- = ap(l - -). 
, we find 



(20) 



a z 5xd x \d) = J\5p) + 0(e 2 



(13) 



p—d 

-l + 2a 2 p{2-p) ~ -1 + C(a 2 ). (21) 

Then substituting the two results (|20| ) and (^lj) into (|lg| ) 
yields the desired effective field theory 

Sx 2 



d t p + ca(l - p)d x p - 



25t 



d 2 p + O(e i ,ea 2 )=0, (22) 
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which is the nonlinear Burgers equation for u = c(p — 1) 
with variable sound speed c s — ca and fixed kinematic 
viscosity v = 

In the special case when a = 0, ( p2] ) reduces to the 
diffusion equation. However, there exists an abnormal 
case when a — 1 where the random bit is fixed to the 
value of 1. The 3-bit model reduces to a 2-bit one that 
does not model the Burgers equation. Its collision term 
is Q = (1 — p + )p- and the equilibrium occupations are 
d + = p and d- = 0. The components of the Jacobian 



matrix are J + 
and J-f — J_ = 



and J_ 
- 1 — p, so (|18[) becomes 



1 + p. Hence d + — cL 



d t p+cd x p- 



Sx 2 (d x p) Sx 2 ( p + 1 



St (1 -p) 2 St \p-l 



d 2 p+0(e 3 ) = 0. 

(23) 



See figure [j] for numerical solutions of (22) and ( |23| ) 



fel28 


fell5; 






fel40S 

\ 




fe640 




1 


fe896 




) 



FIG. 1: The 2-qubit quantum algorithm with 8 = \ and £ = £ = 
(solid curve), 3-bit Boghosian-Levermore classical algorithm with 
a = .707 (dashed curve), and the 2-bit classical algorithm (dotted 
curve) with a = 1. There is good agreement between the quantum 
algorithm and the Boghosian-Levermore algorithm in this case for 

c 2 

equal kinematic viscosities v = -^t- However, the 2-bit classical 
algorithm does not model the Burgers equation demonstrating that 
at least 3-bits per node are required in the classical case. Hence 
the quantum algorithm requires less memory to achieve the same 
result. 

Now we consider the quantum algorithm for the Burg- 
ers equation. Initially we encode the qubits with their 
respective occupation probabilities 
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FIG. 2: The 2-qubit quantum algorithm with 9 = 1.5 radians and 
£ = £ = (solid curve) and the 3-bit Boghosian-Levermore classical 
algorithm with a = 0.5 (dotted curve). This demonstrates that the 
quantum algorithm can model a low viscosity fluid when 9 ~ ^ and 
is therefore computationally efficient. 



ignoring the possibility 
gle. The initial ket 



of an internal 



an- 



phasc 

I \q- ) = 

y/ mr\U) + ^ P+ (l-p_)|10) + y/(l- P+ ) P _\01) + 

— p+){l — p-)\00) is transformed by application of 
a unitary matrix: 



W) = uW). 



The identity matrix and n — 



1 




(25) 



denote the single 



qubit number operator, the multi-qubit number opera- 
tors fbi = h 1 and hi = 1 (g> n are used to determine 
the new probabilities of the respective updated qubits 



V- 



= ^'\n\\i>') =p+ + Cl QLG (p + ,p-) (26) 
= (ip'\n 2 \ip') =p- - fl QIjG (p + ,p-). 



\q±) = V^l^ + v 7 ! - p±|o>, 



(24) 



( Pq ) implicitly determines the functional form of the 
mesoscopic collision term f2 QLG (p + ,p_) associated with 
the microscopic operator U . We use a conservative col- 
lision operator as our 2-qubit quantum gate entangling 
the qubits using only the microscopic states |01) and 1 10) 
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(I 



U 











e iC sin 9 



-e~ iC sin 9 e'^cosd 

\o o o iy 



(27) 



Inserting (|27]) into (|2^) , and then substituting the result- 
ing \ip') into (]26|), we find the collision term f2 QLG is H 



^qlg = sin 2 6>(p_ 



(28) 



+ sin 29 cos (C - y/p+ ( 1 - p+ >_ ( 1 - p_ ) 



Again, we compute the components of J: 



on 



dp± 



And this implies 



Vd+(l-d+)d-(l-d-)y 
(35) 



J+ - J_ = -2 sin 2 9(1 + a 2 /), 



(36) 



The equilibrium condition il QLG 



p—d 



l-d+ 1-d- 



= 2cot#cos(C-£)i 



becomes: 



1 - d+ 1 - d- 



(29) 

which is a legitimate statement of detailed-balance of col- 
lisions at the mesoscopic scale since the quantum model's 
evolution operator is unitary. We take the equilibrium 
occupation probabilities to have the following form: 



d+ = 



1 



and 



1 



jz - 



* + 1 

7 



(30) 



Substituting (|30|) into ( p9j ) gives a quadratic equation in 7 
that has the solution 7 = y/a 2 + 1 + a or i = vc?+l — 
a, where a = cot6>cos(£ — £). Next, substituting ( |30| ) 
into the total number density, p = d + + <i_ , we obtain a 
quadratic equation in z 



/>■- ( 7+ i ] [i, - 1): + /> - 2 = 0. 



(31) 



Substituting the positive root solution of ( pl[ ) into ( |30[ ) 
we find after much algebraic manipulation 



where the factor / = f(a,p) is too complicated an ex- 
pression to write out here but has the important property 
that /(a, p) = 1 + 0(a). Finally, substituting the two 
results ( |34| ) and ( |36| ) into (jl^) gives the effective field 
theory 



dtp+c cot9cos((—{;)(l—p)d x p — cot 2 



^dlp+0(s\ea% 
(37) 

which is the nonlinear Burgers equation for u — 
c(p — 1) with independently tunable sound speed c s — 
ccot6*cos(C — £) an d kinematic viscosity v = cot 2 6$fj- 
by appropriately choosing the Euler angles in fl27|). In 
our quantum case, the trigonometric term cot 9 cos(£ — £) 
plays the role of the expectation value a of the additional 
random bit required in the classical model. Figures [l] and 
H show the time evolution of the 2-qubit quantum algo- 
rithm versus the 3-bit classical Boghosian-Levermore al- 
gorithm both carried out on a N x = 256 lattice. The 
vertical axis is the particle number density p = p + + p_ 
plotted in the range of 0.5 < p < 1.5. The time step is 
in the upper left corner of each snapshot. The viscosity 
of the quantum model is close to zero. 



+ 7 2 + (1 - 7> - 7A/ + 7) (P - I) 2 + 4(p - 2) P< 



The quantum algorithm is unconditionaly stable, 
obeys detailed-balance, requires less memory than its 



2(1 - 7 2 ) 



Then substituting 7 
result 



± 2 ^ 2a 



1 + a into 03 



1 + a 2 - v/l + a 2 (p- l) 2 



This implies that 

d+ - d- = -- (Vl + a 2 - Vl + a 2 ^- l) 2 



classical counterpart, and can achieve arbitrarily high 

(32) Reynolds numbers. Having a variable transport coeffi- 
gives the cient that can be made small, it is consistent with the in- 
viscid Burgers equations when the Euler angle 9 ~ ^. It 
is possible to generalize this type of quantum algorithm 

(33) ^° three-dimensions to efficiently handle the important 
application of computational fluid dynamics. 

I would like to acknowledge Professor Boghosian for 
helpful discussions about his classical lattice-gas model 

(34) and Owen Cote for helpful discussions about turbulence. 
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